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Qh| Abstract. We construct an explicit representation of viscosity solutions of the 

■ Cauchy problem for the Hamilton- Jacobi equation (H, a) on a given domain Q, = 

(0, T) x W 1 . It is known that, if the Hamiltonian H = H(t,p) is not a convex (or 

■ concave) function in p, or H{-,p) may change its sign on (0,T), then the Hopf- 

type formula does not define a viscosity solution on VL. Under some assumptions 
for H(t,p) on the subdomains (i^, ti+i) x K." C £1, we are able to arrange "partial 
solutions" given by the Hopf-type formula to get a viscosity solution on £1. Then 
we study the semiconvexity of the solution as well as its relations to characteristics. 



Oh 
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1. Introduction 



This paper is devoted to constructing a representation formula for viscosity solu- 
00 . tions of the Cauchy problem for the Hamilton- Jacobi equation (H, a) of the form 

5v (i- 1 ) — + H(t,x,D x u) = 0,(t,x)eQ = (0,T)x 

O 

(N| ■ (1.2) u(0,x) = a(x) , x G 



The classical Hopf-Lax-Oleinik formula plays an important role in studying prop- 
erties of solutions of problem (1.1)-(1.2), where H = H(p). Some generalized versions 
of these formulas have been established, and they are useful tools for investigating 
problems of variational calculus, differential games, etc. In most cases, the data 
used to construct the formulas are concerned with the convexity (or concavity) in 
the global setting. This may meet requirements of calculus of variation or opti- 
mal control problems. Until now, no representation formula of global solutions of 
Hamilton- Jacobi equations without concerning the convexity/ concavity or its related 
versions has been constructed. Nevertheless, in the theory of PDEs, the convexity 
posed on the given data is not an obligation in nature. 

In this paper, we consider nonautonomous Hamilton- Jacobi equations in a new 
situation. We suppose that the behavior of the Hamiltonian H(t, x,p) may vary on 
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each subdomain x K" x 1" of the domain [0, T] x R n x IR n , where 



For example, H(t, x, •) may be convex for all (t, x) G [ti-i, U] x R n , and then H(t, x, •) 
changes into a concave function in the successive subdomain (t, x) G x R n . 

The paper is structured as follows. In section 2, we introduce a way to joint 
viscosity solutions of problem (1.1)-(1.2) on subregions of the form [tj,t i+1 ] x R n in 
order to obtain a global solution called a "layered viscosity solution" of the problem 
as a whole. In section 3, we examine the case where the Hamiltonian H = H(t,p) 
is continuous and the initial data <j(x) is convex. Under some assumptions on the 
given data, we use the Hopf-type formula for "partial solutions" 



to build a "layered viscosity solution" of the Cauchy problem on the whole domain. 
We also prove that the obtained viscosity solution inherits the semiconvexity analo- 
gous to "partial solutions" . In section 4, we study relations between the Hopf-type 
formula and the "layered viscosity solution" as well as its singularity based on the 
semiconvexity of the solution u(t,x) in connection with the characteristics. 

Our notations are standard in the field of nonlinear PDEs. For a positive real 
number T, denote Q = (0, T) x R n . Let | . | and (., .) be the Euclidean norm and the 
scalar product in R n , respectively. For a function u : Q — > R, we denote by D x u 
the gradient of u with respect to variable x, i.e., D x u = (u xi , . . . , u Xn ). 

A continuous functions u : [0, T] x R™ — > M is called a viscosity solution of the 
Cauchy problem (1.1)-(1.2) on Q provided that the following hold: 

(i) u(0,x) = a(x) for all x G R n ; 

(ii) For each v G C 1 (fi), if u — v has a local maximum at a point (to, Xq) G Q, then 

v t (to, £o) + H(t , x , D x v(t , x )) < 0, 
and if m — v has a local minimum at a point (to, xq) G Q, then 

v t (to, x ) + H(t , x , D x v(t , x Q )) > 0. 



Consider the partial differential equation F(x,u, Du) = on a domain O C M m . 
Suppose that O can be divided in two open subsets 0\ and O2 by a C l surface T, 
that is O = 0\ UO2 UT, then several compatible conditions can be added on Oi and 
T, as well as C^-solutions Ui of the equation F(x, u, Du) — on Cj U V for i — 1, 2, 
to get a viscosity solution on whole O; see [31 Theorem 1.3]. However, in the case 
of the Cauchy problem, the situation seems to be rather simple as follows. 

Theorem 2.1. Let H(t,x,p) be a continuous function on [0, T] x R™ x R n and let 
o~i(x), o~2{x) be continuous functions on W 1 . For ti G (0,T) ; consider the following 



= t x < t 2 < ■ ■ ■ < t k = T. 



uAt, x) = max {(x, 




2. 



Layered viscosity solution" in general cases 



Cauchy problems: 
(2.1) 
and 
(2.2) 



u t + H(t, x, D x u) = 0, (t, x) G (0, ti) x 
U(0,X) = CTi(x), xel", 



r, 



Mt + H(t,x,D x u) = 0, (t,x) G (*i,T) x 
u{ti,x) = a 2 (x), x G R™. 

Suppose that problems (2.1) and (2.2) /jave viscosity solutions Ux(t, x) and U2(t, x) 
on the domains Qi = (0, tx) x R n and Q 2 — {ti,T) x ~R n , respectively, where cr 2 (a;) = 
x) /or a// x G R™. T/ien i/ie Cauchy problem 

(u t + H(t, x, D x u) = 0, (t, x) G Q = (0, T) x R n , 

( ^ 2 ' 3 ' ) \n(0,x) = 0-1 (x), xGR" 

/ias a viscosity solution of the form 

ux(t,x), (£, x) G [0, ti] x R n , 



(2.4) u(t,x) 



u 2 (*,aO, (t,x) G [ti,T] x 



Proof. Under the assumptions made, the function u defined by (12. 4p is continuous. 
We will use the "extrema at a terminal time" property of the Cauchy problem. Let 
(to,xo) G £1 and v G C 1 (fi). Assume that u — v has a local maximum at (£o,#o)- If 
(£0,^0) ^ * = 1)2, then the restriction of u(t,x) on fij coincides with Ui(t,x). 
Thus, the inequality 

v l t (t , x ) + F(£ , aro, D x v\t , x )) < 0, 

holds since uAt^x) is a viscosity subsolution on flj, and t> 1 = t>L is also a test 
function such that itj — t>* attains a local maximum at (to,xo). It follows that 

v t (to, x ) + H(t , x , D x v(t , xq)) < 0. 

Now, let to — £1 an d let n — f have a local maximum at (£i,xi) G Q. It is obvious 
that u± —v also attains a local maximum at (tx,Xx) on the set (0, tx] x R n . Arguing 
as in the proof of the [31 Lemma 4.1], we see that 

Vt{to, x ) + H(t , x , D x v(t , x )) < 0, 

that is, u(t, x) is a viscosity subsolution of the problem. Similarly, we are able to 
verify that u(t, x) is a viscosity supersolution. □ 

We call a solution of problem (1.1)-(1.2) given by (2.4) a layered viscosity solution 
of the problem. In the lemma below, we will show that a translation of coordinates 
does not change the properties of viscosity solutions. 

Lemma 2.2. Suppose that H(t,x,p) is a continuous function on [£i,T] x R™ x 
R n , u(t,x) is a continuous function on [t^TjxR", ando~(x) is a continuous function 
on R n . Let 

Hx{r,x,p) = H(tx + r,x,p), Ux{t,x) = u(tx + t,x),t G [0,T — tx] 
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for x,p G R n . Then u(t,x) is a viscosity solution of the Cauchy problem (H,a) 
on (ti,T) x R n if and only ifui(r,x) is a viscosity solution of problem (Hi, a) on 
(0,T-*i) xR n . 

Proof. Observe that if v(t,x) G C^Qt^T] x R n ), then Vx(t,x) = + t,x) is of 
class C 1 ([0,T — 1±] x R n ). Thus, the proof of the lemma is straightforward using the 
definition of a viscosity solution. □ 

Remark 2.3. Suppose that the Hamiltonian H = H{t,x,p) : O — > R is continuous, 
where O = [0,T] x R n x M n . Let t x G (0,T) and put d = [0,ti] x M n x M n , 
£> 2 = [^,T] x M n x M n . Let ^ (resp. // 2 ) be the restriction of H on O x (resp. C 2 ). 
If H(t, x,p) is a convex or concave function with respect to the third variable for all 
(t,x) G [0, T] x R n , a viscosity solution of the problem (H,o~) can be obtained by 
methods of variational calculus under some assumptions; see [2] and the references 
therein. If the convexity of H(t, x, p) in p is changed, for example, Hit, x, p) is convex 
on 0\ and concave on C 2 , then the afore- mentioned method cannot be applied on 
the whole domain [0,T] x M n . Nevertheless, we can get a layered viscosity solution 
for problem (H,a) using Theorem 2.1 



3. "Layered viscosity solution" based on Hopf-type formulas 

In this section, we pay attention to the case where the initial data a(x) is convex 
and H = H{t,p) is a continuous function satisfying some additional conditions. 

3.1. Assumptions. We now consider the Cauchy problem for the Hamilton- Jacobi 
equation: 

du 

(3.1) — + H{t, D x u) = , (t,x)en= (0, T) x W 1 , 

(3.2) u(0,x) = a(x) , x eR n , 

where the Hamiltonian H(t,p) is of class C([0,T] x R n ), and a(x) G C(R n ) is a 
convex function. 

Let a* be the Fenchel conjugate of a. By definition, 

a*(x) = sup{(x,q) - a(q)}. 

q£R n 

We denote by 

D = doma* = {y G R n | a*(y) < +00} 

the effective domain of the convex function a*. The Hopf-type formula for the 
Cauchy problem (3.1)-(3.2) is given by 

(3.3) u(t,x) = max{(x,g) - a*(q) - [ H(T,q)dr} = (a* + [ H(r, -)dr)*(x). 

i£ Rn Jo Jo 

To prove that the Hopf-type formula (3.3) is a viscosity solution of problem (3.1)- 
(3.2), we use the "consecutiveness property" (or "C- property" for short) of viscosity 



solution (see [6]): for all t, s such that < s < t < T, 

(a*+ [ H{T,-)dT)**+ [ H(r,-)drY = \a* + [ H{r,-)dr 
Jo Js J L J 

This relation is true if condition (A2) (to be formulated later) for H(t,p) holds. It 
concerns mainly with the assumption that the Hamiltonian H(-,p) does not change 
its sign on (0,T). If this condition is violated, then the Hopf-type formula (3.3) is 
no longer a viscosity solution of problem (3.1)-(3.2); see also [HIE]. Therefore, we 
will use Theorem 2.1 to overcome the difficulty Suppose that H(-,p) may change 
its sign several times on (0,T). Then on each subinterval (ti,t i+ i), where H(-,p) 
does not change its sign, we use Hopf-type formula to get a viscosity solution Ui(t, x). 
Then we arrange consecutively these "partial viscosity solutions" to obtain a layered 
viscosity solution of the problem on whole interval (0, T). This solution also satisfies 
the semiconvexity as in the case of the Hope-type formula. 

Let us recall some results necessary in further. First, we need the following as- 
sumptions on the domain [a, b] x R™, a < b: 

(Al) For every (to, xq) G [a, b] xW 1 , there exist positive constants r and N such that 
(x,p) — <J*{p) — \ H(r,p)dT < max{(x,g) — cr*(<?) — / H(t, q)dr}, 

Ja \l\< N Ja 

whenever (t, x) G [a, b] x IR n , \t — t \ + \x — x | < r and \p\ > N. 
(A2) H(t,p) is admitted as one of two following forms: 

a) H(t,-) is a convex (or concave) function for all t G [a, b]. 

b) H(t,p) = g{t)h(p) + k{t) for some functions g, h, k, where g{t) does not change 
its sign for all t G [a, b]. 

For each (t, x) G £1, let £(t, x) be the set of all p G W 1 at which the function 



q h-> (x, q) - a*(q) - / H(r, q)dr 

J a 



attains a maximum. In virtue of (Al) on [a, b] x M n , one has £(t, x) ^ 0. 
Then we have the following theorem. 

Theorem 3.1. ([6]). Assume that (Al) and (A2) are satisfied on [0, T] xl n . More- 
over, let H(t,p) G C 1 ([0,T] x ]R n ). Then the Hopf-type formula u(t,x) defined by 
( 13.31) is a viscosity solution of the Cauchy problem (3.1)-(3.2). 

As a direct consequence of Lemma 2.2, the following proposition is used to define 
the Hopf-type formula for the initial time to 7^ 0. 

Proposition 3.2. LetO <t\ <T. Suppose thato~(x) is a convex function onMJ 1 , and 
H(t,p) is a continuously differentiable function on [t±,T] x IR n satisfying conditions 
(Al) and (A2) on [t h T] x R n . Then the function 

u(t,x) = (a*+ [ H(r,-)dr)*(x) 
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is a viscosity solution of problem (H,o~) on (ti,T) x R n . 

The following theorem defines a viscosity solution of the Cauchy problem for 
Hamilton- Jacobi equation (if, a) where condition (A2) holds on each subinterval of 



Theorem 3.3. Let g{x) be a convex function on R n and H(t,p) G C 1 ([0,T] x R n ). 
Suppose that there exists t\ G (0, T) such that on each subdomain [0,ti] x IR n and 
[t 1,1*] x R n 7 conditions (Al) and (A2) are satisfied, where g(-) is replaced by u(ti, ■) 
for condition (Al) on [t\,T] x R n . Then the function u(t,x) given by 



is a viscosity solution of the Cauchy problem (H,o~) on (0, T) x R™. 

Proof. We note first that, for fixed t\ G (0, T), the function uu(-) — u(ti, •) is convex. 
By Theorem 3.1 and Proposition 3.2, and then applying Theorem 2.1, we deduce 



Remark 3.4. From Theorem 3.3, we have several corollaries. For example, we can 
suppose that the interval [0, T] can be split into m subintervals [tj,tj + i] on which 

(i) H(t,-) is convex (or concave), or 

(ii) If H(t, x) = g(t)h(p) + k(t) then g(t) does not change its sign on each subinterval. 

Of course, on [0, T], the Hamiltonian H(t,p) can be mixed by forms (i) or (ii) from 
one subinterval to another. Then a viscosity solution u(t,x) on [0, T] x R n can be 
defined as m— "layered viscosity solution" via the solutions Ui(t,x) of the problem 
(H, Gi) on (ti, t i+ i) xl" given by Hopf-type formula, where (Ji(x) = «»(£«, x), g (x) = 
g(x) are convex functions. 

3.2. The semiconvexity of layered viscosity solutions. 

Definition 3.5. A function v : Q — > R is called semiconvex with linear modulus if 
there is a constant C > such that 



for any y\ = (ti,Xi) and y<i = (£2,^2) in ^ and for any A G [0, 1]. 

The function v is semiconcave if and only if — v is semiconvex. In this paper, 
we therefore use some properties of semiconvex functions transferred from those of 
semiconcave functions given in [2]. The theory of semiconcave functions has been 
developed since the last decades of previous century. The reader is referred to the 
aforementioned monograph [2] for a systematic development of the theory. 

In our recent paper [6], we prove that the viscosity solution of the Cauchy problem 
(H, g) given by Hopf type formula (3.3) is a semiconvex function on the domain Q. 
For the "layered viscosity solution" , we also get the same property as in the theorem 



the (0,T). 




n 



the conclusion of the theorem. 



□ 



v(\ yi + (1 - \)y 2 ) - Xv(yi) - (1 - X)v(y 2 ) < A(l - X)-\ Vl - y 2 
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below. In this theorem condition (Al) is automatically satisfied since a is convex 
and Lipschitz on R™, which implies the domain of a* is bounded; see [9]. 

Theorem 3.6. Suppose that H(t,p) G C 1 ([0,T]xR n ) and a is convex and Lipschitz 
on R n . Let t + G (0, T) be a positive real number such that on each subdomain 
[0,t + ] x R™ and [t + ,T] x R™ ; the Hamiltonian H(t,p) satisfies condition (A2). Then 
the layered viscosity solution u(t,x) defined by 



is a semiconvex function on Q. 

Proof. For (t 1 ,Xi) and (t 2 ,x 2 ) i n where t\ < t 2 , and A G [0,1], put (t ,x ) = 
X(t 1 ,x 1 ) + (1-X)(t 2 ,x 2 ). 

Consider the following cases: 

Case 1. < ti < t 2 < t+ or t+ < h < t 2 < T. 

Let D = dom a*. Arguing as in the proof of Theorem 2.4 in [6], we obtain the 
inequality below 

u(t ,x ) < Xu^Xi) + (1 - X)u(t 2 ,x 2 ) + A(l - X)M\t 2 - ti| 2 , 

where M = sup^Q^ p£D \H t {r,p) \ by applying the mean value theorem. Note that 
D is bounded since o~(x) is Lipschitz on R n . 
Case 2. < t x < t < t+ < t 2 < T. 

Take p G £(t , xq) and let 



(3.4) 





(t,x) G [t+,T] x W 



n 



$ = u(t ,x ) - Aw(ti,xi) - (1 - X)u(t 2 ,x 2 ). 



Observe that 




since v** < v for any function v. 
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It follows that 

$ < ( Axi + (1 - X)x 2 ,p) - a*(p) - [ ° H(r,p)dr 

Jo 

- X(x 1 ,p) + Xa*(p) + X H(r,p)dr 

Jo 

- (1 - X)(x 2 ,p) + (1 - X)u*(t + ,p) + (1 - A) / 2 H(r,p)dr 

Jt + 

< A ( jT if (r, p)dr) + (1 - A) ( jf *° if (r, p)dr - jf *" if (r, p)dr) 



+ (l-A)(ix*(f + ,.)-<7*(0)(p) 



< A 



, ( jf 1 H(r, p)dr^j + (1 - A) ( jf 2 if (r, p)dr' 

Therefore, we have for some rj" G [ti, t ], r 2 * e [to ,t 2 ], 
$ < A(t x - t )H(r*,p) + (1 - A)(t 2 - t )H(r;,p) 

< A(*x - (A^ + (1 - A)t 2 ))ff (r^p) + (1 - A)(t 2 - (Atx + (1 - A)t 2 ))ff (r*,p) 

< A (l - A)(ti - t 2 ){H(r*,p) - H{T*,p)) 
<X(l-X)\ tl -t 2 \\H(T*,p)-H(r;,p)\ 

< A(l — A) M |ti - t 2 | 2 , 

where M is defined in Case 1. 
Consequently, 

u(*o, x ) < Au(*i, an) + (1 - X)u(t 2 , x 2 ) + A(l - A)M|ti - t 2 | 2 . 

Case 3. < ti < t+ < t < t 2 < T. 

We can assume that t+ < to because the other possibility has been considered 
in Case 2. Since to = Ati + (1 — A)t 2 , one has A = . The line segment joining 
(ti,xi) and (t 2 ,x 2 ) intersects the plane t — t+ at (t+,x+). Let t = /it + + (1 — /i)t 2 
and t+ = ati + (1 — a)t . Then 

to c\ h — to , t n — 1+ 

(3.5) /i = and a 



t 2 — t+ to — ti 

A simple calculation shows that 

A= afI and 1-A= 1 " " 



1 — fi + afi 1 — /x + a// 

Applying Cases 1 and 2 in the sequel, we have 

«(to, £ ) ^ A* u (*+) + (1 - fj)u(t2, X 2 ) + //(l - yU)M|t 2 - t + | 2 

and 

■u(t + , x + ) < ctu(ti, x\) + (1 — a)u(t , x ) + a(l — a)M\t — t\\ 2 . 



Therefore, 

u(to, xq) < cefiu(ti, Xi)+n(l — a)u(to, Xq) + fj,a(l — a)M\t — *i| 2 

+ (1 - (Jt)u(t 2 , x 2 ) + /i(l - n)M\t 2 - t + \ 2 . 

Thus, 

aji , . 1 — a , . 

«(*o,^o) <i ; + t ; u(t 2 ,x 2 ) 

1 — /i + a/i 1 — /i + a/j, 

+ " Ml = a) M\t - * x | 2 + M\t 2 - * + | 2 

1 — fi + a/J, I — jj + an 

< Xu{ti,xx) + (1 - X)u{t 2 ,x 2 ) + M(A(1 - a)|*o - *i| 2 + n(l - X)\t 2 - t + \ 2 ). 

Let A = A(l — a) |t - *i| 2 + M 1 - A) 1*2 - *+| 2 - Using ([33]) and the fact that 
A = t-t-, we have 

A = A(l - a)(l - A) 2 |* 2 - tx| 2 + — |* 2 - * | 2 

= A(l - A) 2 (l - a) |t 2 - *i| 2 + (1 - A)A-|* 2 - *i| 2 

// 

= A(l - A)|* 2 - *!| 2 ((1 - a)(l - A) + -) < 2A(1 - A)|* 2 - *x| 2 

/i 

A * 2 — *+ r , 

since — = G (J, 1 . 

/x * 2 - *i 

Finally, we get 

u(to, x ) < An(* 1; Xl ) + (1 - A)u(* 2 , x 2 ) + A(l - \)2M\t x - t 2 \ 2 . 

It follows from the definition that the function u(t,x) is semiconvex in Q. The 
theorem has been proved. □ 

Remark 3.7. 1. Under the assumptions of Theorem 3.6, the viscosity solution «(*, x) 
given by H3.4f) is a semiconvex function. Thus, it is also a locally Lipschitz function 
on Q. For * G [* + ,T], we can rewrite u(t,x) as follows: 

u(t, x) = (V + jf + (r, •)*")** + jf F(r, -)<fr 



On the other hand, the Hopf-type formula for the Cauchy problem (3.1)-(3.2) is 
of the following form: 

ft 

(x) 



u H (t,x) = (** + [ H(r,-)dr)*( 
Jo 



The function uu(t,x) is defined on [0,T] x R n , and it is also a locally Lipschitz 
semiconvex function; see [6, 10]. Note that 

{cr*+ f + H(r r )drY* < a* + f + H(r,-)dr. 
Jo Jo 

Therefore, u H (t,x) < u(t,x) for all (t,x) G [t+,T] x W 1 . 
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If the "C-property" for H(t,p) and a holds on [0,T], then u(t,x) = UH(t,x) on 
VL. Otherwise, for some (t, x) G Q, where t > t + , one has u(t, x) > Un(t, x). 

2. It is worth noticing that both functions u(t,x) and un{t,x) are Lipschitz 
solutions of the Cauchy problem (3.1)- (3. 2), (i.e., the solutions are locally Lipschitz 
and satisfy the equation at almost all points of the domain) and moreover, they are 
semiconvex functions. Therefore, the criteria "semiconvexity" or "semiconcavity" 
are not strong enough to get the uniqueness of Lipschitz solution of the Cauchy 
problem for Hamilton - Jacobi equation if the Hamiltonian H(t, x,p) is not concave 
or convex in the variable p; see [4, p. 132]. 

4. Relationship between viscosity solutions and characteristics. 

Singularity 

4.1. The Hopf-type formula and characteristics. Next, we are concerned with 
the function u(x) = u(ti,x) as the initial data for problem (H,u) on the inter- 
val [ti,T]. We will investigate the relationship of characteristic curves between two 
problems {H,a) and (H,u). To this aim, let us recall the Cauchy method of char- 
acteristics for problem (3.1)-(3.2). Note that, to use the method of characteristics, 
at least the initial data is assumed to be of class C 1 . Therefore, some sufficient 
conditions for u(t,x) to be of class C 1 in [0,ti] x M™ will also be presented here. 

From now on, we will use the standing assumptions that H(t,p) G C' 1 ([0, T] x M n ) 
and a G C X (R"). Moreover, a is convex and Lipschitz on R". 

The characteristic differential equations of problem (3.1)-(3.2) are given as follows: 

(4.1) x = H p ; v= (H p ,p)- H ; p = 
with the initial conditions: 

(4.2) x(0) = y ; v(0) = a(y) ; p(0) = a y (y) , y G W\ 

Then a characteristic strip of problem fl3.1|) - fl3.2p (i.e., a solution of the system of 
differential equations (14. ip - (14. 21) ) is defined by 

x = x(t,y) = y+ H p (r,a y (y))dT, 
Jo 

(4.3) < /■* f* 

v = v(t,y) = a(y)+ (H p (T,a y (y)),<7 y (y))dT - H(T,a y (y))dr, 
Jo Jo 

^p = p(t,y) = a y (y). 

The first component of solutions (14.31) is called the characteristic curve (briefly, 
characteristics) emanating from y i.e., the curve defined by 

C : x = x (t,y) = y+ [ H p {r, a y (y))dr, t G [0, T\. 

Jo 

Let (to^o) £ ^- Denoted by £*(to,xo) the set of all y G M. n such that there is 
a characteristic curve emanating from y and passing the point (to,Xo). We have 
£(t ,x ) C a y (£*(t ,x ))] see [8]. Therefore, t(t ,x ) ^ 0. Moreover, 
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Proposition 4.1. ([7]) Let (t ,x ) G fi. Then a characteristic curve (C) passing 
(t , x ) has form 

x = x(t, y) = x + / H p (t, a y (y))dr 

J to 

for some y G £*(t ,x ). 

We say that the characteristic curve C is of type (I) at the point (t ,x ) G f2 if 
<Jy(y) = p G £(to,^o)- If °~y(y) £ °"y(^*(^0) ^o)) \ £(Jo, x o), then C is said to be of type 
(II) at (t ,x ). 

For ti G (0, T), we reproduce in the theorem below a result in [7] stating that the 
viscosity solution u(t,x) of problem (H,cr) is of class C 1 ([0,ti] x R"). 

Theorem 4.2. Let H(t,p) satisfy condition (A2) on [0,*i] x R n and fe* <7 G C^R") 
&e convex and Lipschitz. Suppose that the function 

uj(x) = u(h, x) = (a* + f H(t, -)dTy(x) 

Jo 

is o/C^R"). Tnen the solution 

u (t,x) = (a*+ [ H(r,-)dr)*(x) 
Jo 

is of class C^tMi] x R n ). 

Let (t ,x ) G {t\,T) x R n . Suppose that the characteristic curve C 

x = x + / H{t, a y (y))dr, 
J t 

where y G £*(t ,x ), intersects the plane t — t± at (ti,xi). We have the following 
theorem. 

Theorem 4.3. Let H(t,p) G C^QO.T] x R") and let a G C^R") 6e a convex 
Lipschitz function. Suppose that for some t\ G (0,T), the function 

uj{x) = u{t u x) = (a* + I' H{t, -)dr)*(x) 

Jo 

is of class C 1 (R"). Then every characteristic curve 

C : x = x(t, z) = z + / H p (T,uj z (z))dT, t E [t±,T] 

Jh 

of problem (H, u>) can be extended backward to form a characteristic curve of problem 
(H, a) starting at (0,y), where y = z — J* 1 H p (t, u> z (z))dr. 

Conversely, let C : x — x(t,y) = y + J ' H(t, a y (y))dr be a characteristic curve 
of type (I) of the problem (H,cr) on [0, T] x R™. Then the restriction of C on [ti,T] 
is also a characteristic curve of problem (H,u) on [t\,T] x R n . 
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Proof. Let C\ : x = x\ + f* H p (t, u> z (z))dT, t 1 < t < T, be a characteristic curve 
of the problem (H,w) starting from (ti,Xi). Since uj(x) = u{ti,x) is differentiable 
at (ti,Xi) under the assumptions made, we see that £(ti,Xi) is a singleton and 
q = u> z (z) G £(ti, Xi). Let C : x = Xi + f* H p (r, q)dr, < t < t±, be the characteristic 
curve of type (I) of the problem (H,a) passing (ti,xi). Then the curve C given by 



x = x(t) 



xi+ f H p {r, q)dr, t G [0,ti] 
xi+ / H p (r,oj z (z))dr, t e [ti,T\, 



is exactly the characteristic curve starting from y = x\ — J* 1 H p (t, q)dr of problem 
(H,a) on [0,T] x R n . 
Conversely, if 

C : x = x(t, y) = y + / H(r,a y (y))dT 

Jo 

is a characteristic curve of type (I) of problem (H, a) on [0, T] x W 1 starting from 
y G K™. Then a y (y) = q G xi), where xi = y+ f^ 1 H(t, a y (y))dr, and g = u z (xi). 
Moreover, C can be rewritten as 

x = x(t,y) = xi + / H(r,u y (y))dT, t G [ti,T], 

that is a characteristic curve of the problem (H, cu) on [t±,T] x IR n starting at (ti, xi). 
The theorem has been proved. □ 

4.2. Singularity of viscosity solution. Since a viscosity solution of the Cauchy 
problem (3.1)-(3.2) is a semiconvex function, we can apply the results on the prop- 
agation of singularity of a semiconvex (or semi concave) functions established in [2, 
p.85]. 

We will define the semidifferentials of a continuous function. Let u = u(t,x) : 
Q ->■ R and let (t , x ) G ft. For (h, k) G R x R n , define 

«(t + /i, x + A;) - -u(t , x ) —ph— (q, k) 



Tip, q, h, k) = 



D + u(t , xq) = {(p,g)elx R n | limsup r{p,q,h,k) < 0}, 

(ft,fc)->-(0,0) 

D~u(to, xq) = {(p,g)elxl n liminf r{p,q,h,k) > 0}, 

(ft,fc)^(0,0) 

where p G R, g G R n . 

The set D + -u(t ,x ) (resp. D~u(t , x )) is called the superdifferential (resp. snfr- 
differential) of u(t,x) at (t ,x ). 

In what follows, we denote by D*u(t, x) clxR" the set of all reachable gradients 
of x) at (t, x). That means (p, q) G D*u(t, x) if and only if there exists a sequence 
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(tk,Xk)k C fl \ {(t,x)} such that u(t,x) is differentiable at (tk,Xk) and, 

(tk,x k ) -> (t,x), (u t (t k ,x k ),u x (tk,Xk)) -» (p,q) as k ->■ oo. 

If w(t, x) is semiconvex, then D~u(t, x) = co D*u(t, x), the convex hull of D*u(t, x); 
see [2, p. 54]. 

Suppose that u(t,x) is the Hopf-type formula (3.3) for the problem (3.1)-(3.2). 
Then it is a locally Lipschitz function on Q. 

Let (to ? ^o) £ ^- Define 

n{t ,x ) = {(-H(t ,q),q) | g G £(t ,^o)}- 

Then a relationship between D*u(t , x ) and the set £(t , #o) is given by the following 
theorem. 

Theorem 4.4. ([7]) Assume that condition (A2) is satisfied on [0,T] x IR n . In 
case that H(t,-) is a concave function, we assume in addition that H(t,-) is strictly 
concave for a.e. t in (0,T). Let u(t,x) be the Hopf-type formula for problem (3.1)- 
(3.2). Then for all (to,x ) G tt, we have 

D*u(t ,x ) = H(t ,x ). 

We also denote XX M ) the set of all singular points of a function u(t,x) (i.e., 
the points at which u(t,x) is not differentiable). Using an important result on the 
propagation of singularities of a semiconcave function presented in [2], we have the 
following: 

Theorem 4.5. Let [0,T] = [0, £i]U[ii, T] and suppose that condition (A2) is satisfied 
on each subdomain [0,ti] x M. n and [t±,T] x R". Let (t ,x ) G XX M ) an( ^ assume that 
the 0— level set of the function g{p,q) = p + H(to,q) does not contain any line 
segment [(y>i, (y>2, 92)] C M n+1 ; (pi,?i) 7^ ($2^2)- Then there exists a Lipschitz 
arc x : [0, p\ — > W l+1 with x(0) = (t ,x ) and x(s) G XX M ) f or an V s e (0>P) an( ^ 
x(s) 7^ (to 5^0) / or s > small enough. 

Proof. Note that the set £(t ,x ) is compact. Thus, the set 

D~u(t ,x ) = coD*u(t ,x ) = coH(t ,x ) 

is also compact. It follows that there exists (po,qo) G D~u(t ,x ) such that 

g(Po,qo) = , max g(p,q). 

Since the 0— level set of function g(p,q) does not contain any line segment and 
D~u(t ,x ) is not a singleton as (t ,x ) G XX M )> 

max g(p 7 q)= max (p + //(to, g)) = Po + H(t , q ) = a > 0. 

(p,q)eD u(t ,x ) (p,q)eD-u(t ,x ) 

Thus, (po,qo) G D~u(to,x ) \ D*u(t ,x ). For any 11 G N, we see that 

l/n + p + H(t ,q ) > a. 
Therefore, (po + l/n,qo) ^ D~«(to,^o)- This means that 

(Po,<7o) G dD~u(t ,x ) \D*u(t ,x ). 
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Applying [2, Theorem 4.2.2], we obtain the desired result. □ 

Remark 4.6. 1. From Theorem 4.5, we do not know whether the singular point 
propagates forward or backward in time. 

2. If we use method of characteristics or Hopf-type formula to define generalized 
solutions of the Cauchy problem of Hamilton- Jacobi equation, we see that the singu- 
lar points of solutions may propagate forward in time t. Nevertheless, if the structure 
of the Hamiltonian (convexity, for example) changes in time, then the singularities 
of solution may cease. 

Example. Consider the following Cauchy problem: 

u t + (2t - l)u 2 x = 0, (t, x) G (0, 2) x R, 
u(0, x) = \x\, x G R. 

Since the Hamiltonian H(t,p) = (2t — l)p 2 , as a function of p, is concave if 
t G [0, 1/2] and convex if t G [1/2, 2], the Hopf-type formula of the problem 

u(t, x) = max {xp — (t 2 — t)p 2 } 
pe[-i,i] 

is not a viscosity solution of the problem; see [7] . 

On the other hand, applying Theorem 3.3, we obtain a viscosity solution of the 
problem as follows: 



u(t, x) 



max {xp - {t 2 - t)p 2 }, (t, x) G [0, ~] x R, 
pe[-i,l] 2 

max {xp - {t 2 - t)p 2 } + -, {t, x) G [-, 2] x 
pe[— 1,1] 4 L 



A simple calculation gives us 



u(t, x) 



( 'x|-(t-i) 2 + i (t, x) G [0, ~] x R, 

X 1 i , I . o 1 

+ | x |< 2 (t--) 2 , -<t<2, 



4(t- |) 2 4' 1 1 - v 2 y ' 2 

1^1 - (* - n) 2 + T) kl > 2(t - ^) 2 , 7T < * < 2. 



2 7 4 7 

At the points we have = |x| + |, x G R. 

It is obvious that (|,0) is a singular point of the solution and ^(§,0) = { — 1, 1}. 
Therefore, D*u{\, 0) = {(0, -1), (0, 1)} and D~u(\, 0) = {0} x [-1, 1]. Then 

dD~u{\, 0) \ D*u{\, 0) = {0} x (-1, 1) ± 0. 

By Theorem 4.5, the singularity of u(t, x), say, at (^, 0) propagates. Nevertheless, for 
t G (o> 2], the Hamiltonian is strictly convex. Thus, the solution of the problem is of 
C l {{\, 2] x R). This proves that the singularity of the solution propagates backward. 
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